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Abstract 

In the reversed field pinch RFX-mod strong electron temperature gradients develop when the 
Single-Helical-Axis regime is achieved. Gyrokinetic calculations show that in the region of the 
strong temperature gradients microtearing instabilities are the dominant turbulent mechanism 
acting on the ion Larmor radius scale. The quasi-linear evaluation of the electron thermal conduc- 
tivity is in good agreement with the experimental estimates. 
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Transport barriers in magnetically confined plasmas represent one of the topical subjects 
of fusion science, whose ambitious aim is the controlled release of thermonuclear energy in 
future power plants. Barriers i.e., regions where the transport of particles and energy is 
substantially reduced, are phenomena characterized by complexity, multi-scale interactions 
and self-organization processes. Understanding them is a challenge for both theory and 
modelling, also due to the variety of magnetic confinement devices - tokamaks, stellarators 
and now reversed-field pinches (RFPs) - where they are encountered. In this Letter we 
identify, for the first time, a key mechanism driving the physics of transport in the barriers 
observed in the RFP device RFX-mod [2]. When the self-organized Single- Helical- Axis 
(SHAx) regimes are achieved js], the plasma approaches equilibria with conserved helical 
magnetic surfaces, normally endowed with internal transport barriers (ITBs) which bind 
the hot core Q|. This new experimental condition raises the question whether small-scale 
transport mechanisms, until now overshadowed by the large-scale MHD turbulence, may 
rule the fundamental properties of the RFP plasma. 

The strong gradients typical of ITBs are reservoirs of large amounts of free energy avail- 
able to trigger microinstabilities. Their effect is to increase radial transport and damp a 
further increase of the gradients. In tokamaks, this has been a main topic of research for 
several decades. Conversely, studies applied to the RFP configuration started only recently, 
mainly focussed on electrostatic Ion Temperature Gradient (ITG) modes. In this respect, 
various approaches agree that present-day RFPs are subcritical to such instabilities, due to 
the strong Landau damping of the modes 0, 

An instability that attracted the interest of the tokamak community in the '70s is the 
MicroTearing (MT) mode 7|. The MT mode is a high-wavenumber drift-tearing mode 
driven unstable by the electron temperature gradient, unlike its long- wavelength counterpart 
which is essentially current-driven. The existence of MT modes leads to chains of magnetic 
islands, whose overlapping brings to local stochastization of magnetic field lines. Hence, 
it was speculated that MT modes might provide an effective contribution to the thermal 
diffusivity through electron parallel motion along the wandering field lines. However, these 
modes are expected to become stable at relatively low collisionality, !/<w*, with v electron 
collision frequency and u* electron diamagnetic frequency [8(. These conditions are usually 
met in current high-temperature devices, which are effectively collisionless in the core. For 
this reason, in tokamaks, MT modes are nowadays studied either in connection with edge 



2 



transport 



9], or in the core of medium-temperature devices 10j. Quite recently, however, 



some broader interest in core MT modes has revived in connection with strong ITBs and 
their possible coupling with other microinstabilities, like the Electron Temperature Gradient 
(ETG) modes ll|. The ohmically-heated RFX-mod presently features a moderate peak 
temperature (1 to 1.5 keV) which, coupled to the strong internal temperature gradients, 
should represent an optimal environment for these modes to grow. Given the long-established 
importance of tearing instabilities in the RFP, we found it natural to consider the connection 
between MT modes and ITBs in this machine. 



We present numerical investigations obtained with the gyrokinetic code GS2 |l2j, modi- 
fied to include the RFP geometry 6|. Before discussing the results, we briefly introduce some 
technical aspects related to code and experiment. Simulations are linear, with fluctuations 
of the magnetic vector potential in both parallel and perpendicular components included. 
MT modes are characterized by a very elongated structure along the field lines; in order to 
generate spatially resolved eigenmodes, the simulation domain extends to 9 G [— 80 n, 80 tt], 
where the poloidal angle 9 plays the role of coordinate along the field line in the ballooning 
representation. The radial domain r/a < 0.8 (r G [0, a], a torus minor radius) is considered 
here, as we are interested in the ITB region. The geometry is assumed axisymmetric for 
simplicity, as described in The analysis is based on the profiles of two reference exper- 
imental cases selected from the database of SHAx discharges, i.e., shots 23977 and 24932, 
whose electron temperatures are shown in Fig. [TJ The electron-to-ion temperature ratio 
is estimated to be T e /Ti ~ 1.5, based on spectroscopic Doppler measurements of partially 
ionized impurities J4]. Accordingly, the total 0, ratio of thermal to magnetic pressure, which 
rules the electromagnetic effects in the gyrokinetic equation, is approximately 1.5% in the 
region corresponding to the ITB for both reference discharges. The collision frequency and 
the electron diamagnetic frequency stay roughly in the range v ~ 10 w* ~ 10 5 s -1 . 

Our study shows that the dominant core instability is of microtearing type, Fig. [TJ Es- 
sentially driven by the higher electron temperature gradient, the discharge 23977 turns out 
to be more prone to this type of instability. The tearing nature of such instabilities is clear 
in Fig. 121 where the typical parity of the mode, odd for the electrostatic potential <fi and 
even for the parallel magnetic vector potential An, is shown. In Fig. [3] mode frequency and 
growth rate are plotted as a function of the perpendicular wavenumber k y multiplied by the 
thermal ion Larmor radius pi, for different radial positions along the profile of Fig. HJa). MT 
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modes are characterized by a propagation in the electron diamagnetic direction (negative 
real frequency in the plot). The second branch in the spectrum for r/a = 0.55 at k y pi > 0.25 
corresponds to modes with even parity in 0, characterized by a very elongated envelope of 



13|, 



the eigenfunction in the parallel direction. Such modes resemble those described in 
essentially determined by the passing-electron response. 

To get a more complete picture of MT mode stability in the RFP, we performed a scan over 
a larger range of physical parameters. The MT mode growth rate is mainly affected by elec- 
tron temperature gradient, collision frequency and total /?. The reference values are assumed 
to be those of shot 23977 at r/a = 0.5, see Fig. (Ha). All quantities are normalized to these 
experimental values, and the single parameters of interest are individually modified keeping 
the others fixed. The normalized logarithmic temperature gradient a/L^ e = —adlogT e /dr 
is varied in the interval [2,5], the total (3 and the collision frequency v between one half 
and twice their respective experimental values. In Fig. Ufa) the growth rate is shown as 
a function of k y pi, for two values of cl/Lt b - The peak of the spectrum is moving towards 
larger k y pi as the gradient is increased, together with the maximum of the spectrum. Frame 
(d) summarizes the dependence of the maximum of the spectrum on a/LT e , showing a good 
linear trend and a stability threshold for a/L^ e ~ 2. Fig. E(b) deals with the f3 parameter 
dependence. Large /3 increases the growth rate and the wavenumber corresponding to the 
maximum growth rate. As /3 decreases, MT modes tend to vanish, and a sub-dominant 
branch, peaked around k y pi ~ 0.3 and not visible for the reference experimental /3, is dis- 
closed; again, such modes possess very extended tails along 6 and qualitatively replicate 
those of 13| . The maximum 7 versus shows a strong increasing trend, see frame (e). In 
view of the development of plasma density control and of the capability to achieve higher j3, 
MT-driven turbulence could become a concern for ITB degradation. In Fig. HJ^c) the spec- 
trum is reported for the two selected values of the collision frequency, with the maximum 
growth rate as a function of v drawn in frame (f). As reported in 8|, the growth rate of 
MT modes is expected to have a maximum for u/u* of the order of tens (depending on the 
value of r] e = d\ogT e /d\ogn, with n plasma density), and to decay logarithmically for larger 
frequencies, until stability is found. Similarly, a few additional cases more distant from 
the scan centered on the reference collisionality of Fig. H(c,f), indicate that ^(v) features a 
maximum just around the experimental value of v/oj* ~ 10. 

We have established that MT modes are unstable over a significant range of wavenum- 
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bers, but substantial heat transport ensues only in the presence of appreciable field line 
stochastization. For this to happen, the Chirikov island overlap criterion adapted to the 
MT scenario 14| yields that unstable modes must have poloidal numbers m greater than 



the threshold m = g/(2|g'|pj) 1 / 2 , with q safety factor and q' its radial derivative. This 
makes the minimum separation between mode rational surfaces, Ar ~ Q 2 /fn 2 \(l'\, less than 
a couple of ion Larmor radii, which is the typical width of the resistive layer. In our case 
mo ~ 3 in the region where MT modes are unstable, whereas the peaks of the mid-radius 
spectra in Fig. [3^b) approximately correspond torn ~ 10, making the inequality m > itlq 
fulfilled for a wide subset of wavenumbers. 

To quantify the heat transport deriving from the saturated state of MT turbulence, the 
electron thermal conductivity \ must be estimated. Because of the difficulty involved in 
calculating such a quantity from direct nonlinear numerical simulations - to our knowledge 
this has never actually been done for this kind of turbulence - people commonly resort to the 
quasilinear estimate by Rechester and Rosenbluth 15] . An upper limit for the conductivity 
is given by the collisionless estimate x — D m v t ^ e , where D m is the magnetic diffusion 
coefficient, ratio of the averaged squared radial displacement to the distance traveled along 
the field lines, D m = ((Ar) 2 )/2Z, and w t h, e is the electron thermal speed. The diffusion 
coefficient itself is usually written as D m ~ b 2 L c , with b = 5B/B normalized magnetic 
field perturbation and L c longitudinal correlation length of the field. The saturated level of 
the magnetic field perturbation is provided by the quasi-linear estimate by Drake et al jsj]: 
b ~ p P lLr„. In tokamaks L c is usually assumed approximately equal to the connection 



length 16(. Due to the low value of q, for an RFP the connection length is expected to be 



much smaller, close to 27ra as the reversal is approached. To get a more accurate evaluation 
of L c we relied on numerical simulations with the field-line-tracing code Nemato [l7]. The 
overall magnetic field is built by superposing a large population of modes with poloidal 
and toroidal wave-numbers (m, n) determined by the spectra of Fig. [31 and with amplitude 
such that their cumulative strength is locally ~ p e /LT e - The resulting magnetic field cross- 
section is shown in Fig. [5]^a). The code allows us to compute the Eulerian spatial correlation 
C{1) = (b(l) b(0)) I (6 2 (0)), averaged over an ensemble of trajectories, as sketched in Fig.[5](b). 
The correlation length, given by L c = f£°C(l)dl, becomes constant for large values of the 
upper limit of integration, where C(l) tightly oscillates around 0. It is confirmed that 
L c 2ira, this value turning out to be slightly decreasing for larger amplitudes of the 
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magnetic field perturbation. In frame (c) we show the direct calculation of ((Ar) 2 ) as a 
function of /. The inferred numerical value of D m agrees with the quantity b 2 L c . When the 
values of L c , or directly the value of D m , are inserted into the expression for the quasilinear 
thermal conductivity, we obtain \ ~ 5 4 20 m 2 /s. This turns out to be compatible with 
the estimate for \ given at the barrier location by power balance analysis j^l, [l^, where one 
recovers figures in the range 5 4-50 m 2 /s. As a side remark, our evaluation of x concerns the 
ITB MT-dominated region only, and does not take into account any possible mixing with 
the residual MHD turbulence acting in the edge of the plasma. 

As a further support to MT turbulence as a major drive of transport across ITBs, we note 
that marginal stability of MT modes occurs at a/L Te > 2, the precise figure depending upon 
a combination of several parameters [Fig. 11(d)] . Consistently with the two sample cases 
addressed in this Letter, first statistical investigations on the RFX-mod database suggest 
that the experimental electron Lx e at the barrier typically exceeds the value a/4 [19J. A 
lower bound of Lj- e is likely to be related to a gradient-driven turbulence mechanism which 
i) starts operating above a given threshold in the gradient a/Lj^; ii) activates turbulent 
structures which increase local transport; iii) forbids reaching higher temperature gradients 
(hence the upper bound of a/Lj- e ), due to the high cross- field thermal conductivity which 
tends to smooth the T e profile itself. 

As a final remark, we notice that our findings are consistent with the results from spherical 
tokamaks lOj, suggesting that a similar physics is running in the background. This is not 



at all surprising, as RFPs are acquiring a level of magnetic order until now deemed peculiar 
to tokamaks or stellarators. 

To conclude, this work suggests that MT modes may contribute significantly to heat 
transport across ITBs in present-day RFX-mod plasmas. Our results are in good agreement 
with the experimental estimates of the thermal conductivity. It is still a matter of inves- 
tigation how much MT modes may affect future RFP performances, with special attention 
paid to the scalings in (3 and collisionality. 
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FIG. 1: (Color online) Electron temperature profiles (square symbols) for shot 23977 (a) and 
24932 (b) at 64 ms, with the growth rate of the most unstable MT mode (circles); the gray-tone 
corresponds to the mode growth rate. The edge region r/a > 0.8 is not considered in the present 
analysis. 
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FIG. 2: (Color online) Normalized fluctuations of electrostatic potential (a) and parallel magnetic 
vector potential (b) for the profile of Fig. QJa) at normalized radius r/a = 0.55 and k y pi = 0.18. 
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FIG. 3: (Color online) Mode frequency (a) and growth rate (b) of the most unstable mode as a 
function of the wavenumber k y pi, for different mid-radius positions along the profile of Fig. [T^a). 
O symbols represent modes which do not have microtearing parity. Frequencies are normalized to 
Vt\i/ a i with u t h,i = (Ti/rrii) 1 / 2 ion thermal speed and a minor radius. 
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FIG. 4: (Color online) Dependence of the growth rate of the most unstable mode on k y pi for 
different values of a/Lr e (a), (3//3 exp (b), and ^/f ex p ( c )' * ne experimental parameters referring 
to the profile shown in Fig. QJa) at r/a = 0.5; O symbols represent modes which do not have 
microtearing parity. Frames (d), (e) and (f) summarize the dependence of the maximum of the 
above spectra versus the respective parameter, with the experimental conditions marked by square 
symbols. 
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FIG. 5: Magnetic field line section at 9 = n, in frame (a), with corresponding longitudinal cor- 
relation function of the normalized radial field perturbation (b), and averaged squared radial dis- 
placement of the field lines as a function of the covered distance (c). Calculations of (b) and (c) 
are based on 2500 field line trajectories. 
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